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ABSTRACT 

Gamma-ray bursts are now known to be a cosmological population of objects, which 
are often accompanied by X-ray and optical afterglows. The total energy emitted in 
the afterglow can be similar to the energy radiated in the gamma-ray burst itself. 

If a galaxy containing a large column density of dust is near the line of sight to a 
gamma-ray burst, small-angle scattering of the X-rays due to diffraction by the dust 
grains will give rise to an X-ray echo of the afterglow. A measurement of the angular 
size of the echo at a certain time after the afterglow is observed yields a combination 
of the angular diameter distances to the scattering galaxy and the gamma-ray burst 
that can be used to constrain cosmological models in the same way as a time delay in 
a gravitational lens. The scattering galaxy will generally cause gravitational leasing as 
well, and this should modify the shape of the X-ray echo from a circular ring. 

The main difficulty in detecting this phenomenon is the very low flux expected for 
the echo. The flux can be increased when the gamma-ray burst is highly magnified 
by gravitational leasing, or when the deflecting galaxy is at low redshift. X-ray echos 
of continuous (but variable) sources, such as quasars, may also be detectable with 
high-resolution instruments and would allow similar measurements. 


Subject headings: dust - gamma rays: bursts - gravitational leasing - X-rays: general 


1. Introduction 

The presence of dust grains along the line of sight to an X-ray point source will result in 
a halo of scattered X-rays around the source, which will be time delayed relative to the X-rays 
that were not scattered. The effect is caused by small-angle scattering of the X-rays by the dust 
grains, due to the diffraction of the wavefront when some of the X-rays are absorbed by the 
dust grain. The process has been discussed previously by several authors (Overbeck 1965; Slysh 
1969; Hayakawa 1970; Triimper Sz Schonfelder 1973; Alcock & Hatchett 1978). The applications 
that were considered were generally restricted to galactic sources. The first observations of the 
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phenomenon were made by the Einstein observatory (Mauche and Gorenstein 1986 and references 
therein). 

In the case of an extragalactic X-ray source with a large column density of dust on the line of 
sight, the dust is likely to be concentrated on an intervening galaxy. Therefore, all the scattering 
takes place on a single screen, implying that the time delay is uniquely determined by the angular 
separation from the point source. Thus, a variable X-ray quasar located behing a dusty galaxy 
should show a variable intensity pattern in the surrounding X-ray halo, with the radial dependence 
reflecting the past lightcurve of the quasar. If the time delay at a given angular separation from 
the source can be measured, a combination of the angular diameter distances to the scattering 
galaxy and the X-ray source is obtained, using a method analogous to the measurement of angular 
diameter distances from time delays in gravitational lensing (Refsdal 1966). This method of 
measuring distances was discussed by Triimper &: Schonfelder for galactic sources; in this case, the 
dust is likely to be distributed along the line of sight, and the past lightcurve of the X-ray source 
should be smeared out in the X-ray halo depending on the dust distribution. 

Recently, gamma-ray bursts have been proved to be of cosmological origin, and the associated 
X-ray and optical afterglows have been discovered (van Paradijs et ah 1997; Djorgovski et ah 
1997; Metzger et al. 1997; Costa et al. 1997). The presence of the afterglow was predicted from 
the generic model of interaction of a relativistic fireball with a surrounding medium (Vietri 1997; 
Meszaros &: Rees 1997; Wijers, Rees, &: Meszaros 1997). Gamma-ray bursts may prove to be a 
useful source for doing this type of observation. 

The paper is organized as follows. In §2 we derive the general characteristics of the scattered 
X-rays: the angular size of the halo, the time delay, and the intensity. In general, gravitational 
lensing effects are likely to be present, since a dusty galaxy on the line of sight will typically cause 
a gravitational deflection of the same order as the scattering angle of the X-rays; this is analyzed 
in S3. Our results are discussed in S4. 


2. Characteristics of the X-ray Echo 

We consider an X-ray source from which, in the absence of scattering, a total fluence (or 
energy per unit area) cq is received instataneously; i.e., the flux from the source can be written 
as F = eQS{t — to), where to is the time when the non-scattered X-rays arrive. This may be a 
good approximation for a gamma-ray burst afterglow. The source is at redshift and angular 
diameter distance Dg, and the dusty galaxy is at redshift Zd and angular diameter distance Dd- 
Any scattered X-rays received at an angle 6 from the source will arrive with a time delay td given 
by: 

, _ + Zd) DdDg 


( 1 ) 
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where D^s is the angular diameter distance from the deflector to the source. This equation is 
generally valid for any cosmological model, and is derived in the same way as the equation for the 
time delay in gravitational lensing (Refsdal 1966). 

Thus, the “echo” of an instantaneous pulse of X-rays consists of a ring of angular size 9 at 
time td after the pulse has been received. 


2.1. Intensity of the X-ray Echo 


To calculate the intensity of the X-ray echo, we need to know the scattering cross section as 
a function of the scattering angle a. We consider spherical dust grains of radius a. The typical 
scattering angle of X-rays of wavelength A is a ~ A/(27ra); therefore. X-rays will generally be 
scattered at very small angles by grains with sizes a ~ l^m. The total scattering cross section 
of a dust grain is given by ax = Q^vra^, where Qs is the scattering efficiency. For the case of 
magnesium silicate, a typical material in dust grains, the scattering efficiency can be approximated 
as (see van de Hulst 1957, Alcock k, Hatchett 1978) 


Qs ^ 
Qs ^ 


0.7 

1.5 



otherwise. 


if Qs < 1 ; 


( 2 ) 


The transition between the two regimes occurs when the phase shift across the dust grain, 

IpI = (47ra)/A \m — 1| (where m is the refractive index of the dust), reaches unity. For magnesium 
silicate, \p\ ~ (A/2 A) (a/1 pm). The detailed form of Qs is shown in Figure 1 of Alcock & Hatchett 
(the numbers in the horizontal top and bottom axes in this figure are erroneous; see Figs. 32 and 
33 in van de Hulst 1957 for similar curves). 


The differential cross section can be written as 


da 

= (JX ■ 

da 


2Txa 


CT T. 


2'Ka 


a 


( 3 ) 


The dimensionless function an is normalized to dxanix) = 1. This expression is valid in 
the limit of small scattering angles (this is why the upper limit in the integral normalizing an 
is replaced by infinity). The function an can be calculated for spherical grains under different 
approximations, depending on the value of |/3|. The result is shown in Figure 2 of Alcock & 
Hatchett; generally, <t„(x) rises linearly with x for x <C 1 (as required by continuity of the 
scattering cross section per unit solid angle for zero deflection), peaks near x ~ 1.5, and then 
declines rapidly and may have oscillations of decreasing amplitude at large x (analogous to the 
familiar diffraction pattern of a circular aperture). Practically all the X-rays are scattered at 
angles with x < 4. Of course, dust grains should not be exactly spherical in practice, and there 
should be a wide distribution of grain sizes. 
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Considering now a variable X-ray source with an unabsorbed flux as a function of time FQ{t), 
and assuming that the fraction of X-rays scattered in the intervening galaxy is small (so that 
multiple scatterings can be neglected), the surface brightness S of the X-ray echo at angular 
separation 9 and position angle (p from the source is 

S{9,4>]t) sinO d6 dip = Fo{t — td) T{9,(j)) — da {d4>/2TT) . (4) 

This equation expresses the fact that the flux received along a direction within (d6, dip) is equal to 
the total flux of the source, times the probability of scattering (equal to t(0, </>), which is assumed 
to be small), times the probability that the photon is scattered in the direction within (da, dip). 
Expressing the scattering angle in terms of the angular separation, a = 9 (Dg/Dds), we have 


S{9,(p-,t) 


Fo{t-td) 
sin 9 


t{9,(P) 


aPs 

^Dds 




2'Ka Dg \ 

\Dds ) ■ 
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When the variable source emits most of its energy on a time much shorter than td (i-e., the flux is 
approximated as To(t) = eo d{t — to) as before), the echo is a narrow ring with total flux F^, given 
by 


Fe{t) 


eo aPs 

td 2XPds 


9 an 


2'KaPg 

^Pds 


r2TT 


d(pT{9, ip) 


( 6 ) 


where the angle 9 is to be calculated from the time delay td 


t — to from equation (||) . 


The optical depth for X-ray scattering is due to dust grains of different sizes, and can be 
expressed as r = / Nd{a) aria) da, where Nd{a) da is the column density of dust grains of radius 
a. The dust column density can be written as Nd{a) = Ni, Zd{a) (Aira^pd/3mp)~^, where Ni, is the 
column density of baryons in the interstellar medium of the intervening galaxy, Zd{a) da is the 
fraction of mass of the interstellar matter in the form of grains of radius a (i.e., the “metallicity” 
in dust grains of radius a), pd is the dust density, and nip is the proton mass. We assume that 
Zd{a) is independent of position. Substituting into equation (^), we obtain 

F,{t) = ^f^^^^9 dcPN,{9,cP) j daZd{a)Qg{\,a) an{^^^9^ . (7) 


In order to proceed further, we will need to discuss the model for the function Zd{a). Studies 
of interstellar reddening indicate that Zd{a) oc for small grains, a ^ 0.1 pm, (Mathis, Rumpl, 
& Nordsieck 1977). The function must turn over at larger size so that the total mass in grains 
converges. While most of the mass seems to be in grains of size near 0.1 pm (Martin 1978; Savage 
& Mathis 1979; Mitsuda et ah 1990), it seems natural that the mass is distributed over a wide 
range of a. Thus, as an example, we shall assume Zd{a) oc for a > Q.lpm. The scattering 
of X-rays in the regime we shall discuss is dominated by grains with a ^ 1 pm, so our estimate of 
the flux is sensitive to the uncertain form of Zd{a), which may vary depending on the interstellar 
environment of the galaxy. For the adopted model, if Zt is the total metallicity in dust grains, we 
have aZd{a) = Zt/4{a/0.1 pm)~^'^ . A chemically evolved galaxy can be reasonably assumed to 
have Zt Zq. 
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Let us now examine carefully the terms appearing in the integral over the dust grain radius 
in equation]^ The function an is sharply peaked at the radius Op given by 2'KapDs6/{XD^s) — 1-5 
(see Fig. 2 of Alcock & Hatchett), and is normalized to / an{x)dx = 1. As long as 
a > 1 ^m (3 A/A) = oi, we can use Qs — 1.5 (see the discussion above; this value for Qs is an 
average over the range 1 ^ a/ai ^ 10, which will be the range of interest). Smaller grains are not 
effective for X-ray scattering because Qs is very small. For the assumed shape Zd{a) oc and 

if Up > ai, a sufficient approximation is obtained by substituting Xd(a) in equation ^ by its value 
at the peak of the function a^ Thus, when Op > oi, the integral can be approximated as 


d(x Zfiici) Qsi^X^ (i) an 


2-KaDs 
^Dds J 


— Qip Zd{cLp') 


Ztl4: 


ap -0.5 


0.1 ^m 


The condition Op > oi implies 
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Ur, ~ 1.5 
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or 

For shorter wavelengths (when the angular separation 0, and therefore the time delay, is fixed), 
the integral rapidly becomes much smaller; essentially, the intensity of the scattered rays is greatly 
reduced because the grains that are small enough to cause scattering at an angle 9 are mostly 
transparent to the radiation. 


Thus, when equation (IT) is satisfied, the flux of the X-ray echo is 


eo dnmpDse - ^ ^ 

Fp — „ , _ N^(^9j dp ZfiyQpj 


~ 0.03 


^d ^Pd^Dds 

eo 1A 0 


Nb 


Z, 


ap -0.5 


td Dds A 1" 1022 cm-2 0.02 0.1/im 


( 11 ) 

( 12 ) 


where Ni,{6) is the mean gas column density over the ring of angular size 0, and we have used the 
value pd = 3gcm“^ (appropriate for magnesium silicate). Notice that A is the wavelength in the 
frame of the deflector in all our equations (i.e., the observed wavelength is A(1 -|- z^)). 

As an example, we consider a gamma-ray burst at redshift 2 ; ~ 1 similar to GRB970228. The 
total X-ray fluence of this burst in the energy band from 2 to 10 keV was 4 x 10“® ergcm“^, about 
40% of the fluence in the gamma-ray burst itself in the band 40 — 700 keV (Costa et al. 1997; van 
Paradijs et al. 1997; Wijers, Rees, & Meszaros 1997). A large fraction of the X-ray fluence was 
emitted in the beginning of the afterglow, about 50 seconds after the burst, and the flux decayed as 
^-1.33 (Qosta et al. 1997). The X-rays that can be scattered most effectively are in the range 2 to 
6 keV (at lower energies, the scattering is large and the X-rays arrive with a very long time delay, 
and at higher energies the condition Op > ai is not obeyed). We take the fluence in this narrower 
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energy band to be 2 x 10 ® erg cm For a time delay of one year, eo/td = 10 ergs ^ cm 


If 


the deflector is at (1 + Zd)Dii ~ 10^ Mpc, and taking Ds/D^s — 2, then equation (|^) gives 9 ~ 3", 
and equations (^) and (0) give A > 2A, Op > 1.5/rm. Using the fiducial values of N^, and Zt in 
equation (12), and A = 2A, we find a flux Fg ~ 10“^^ ergcm“^ s“^. At Zg = 1, the corresponding 
luminosity is Lg = 47r(l + Zs)'^DgFg ~ 0.02Fio/td ~ 5 x lO^^ergs”^. This luminosity is much larger 
than that of normal galaxies, so the X-ray echo can be much brighter than any intrinsic emission 
from the intervening galaxy. 


We see that the dust grains that can produce these X-ray echos are generally of a large 
size compared to most of the grains, which are responsible for interstellar reddening. Small dust 
grains should of course give rise to X-ray echos of larger angular size, but these are more difficult 
to observe because of the longer time delay they imply, and because most random intervening 
galaxies are not larger than a few arc seconds. 


If the X-ray echo can be resolved, and the angular radius of the ring is measured at a time 
after the known occurrence of the gamma-ray burst, then equation (|^ can be used to calculate the 
value DfiDs/Dds- Assuming that the redshifts Zd and Zg are known, this gives a useful constraint 
on the Hubble constant, as well as the value of the radius of curvature and cosmological constant 
of the universe. These can be constrained separately once the observation is done in several 
sources at different redshifts. If the X-ray echo is not resolved, one might still be able to infer 
the angular size of the ring by studying the intervening galaxy in optical and infrared light, and 
inferring the distribution of dust. The lightcurve of the unresolved X-ray echo might then reveal 
the region of the galaxy that the echo is moving through. For example, if the intervening galaxy 
has a smooth dust density distribution peaked at the center, but the position of the gamma-ray 
burst is displaced from the center of the galaxy, then the time when the X-ray echo peaks in 
intensity should correspond to the time when the echo has passed through the center of the galaxy. 


The main difficulty in detecting the X-ray echo will be the very low flux expected, 

~ 10“^® ergcm“^ s“^. For example, the few X-ray clusters of galaxies detected so far at high 
redshift have fluxes above 10“^^ erg cm“^ s“^ (Luppino &: Gioia 1995) (however, the X-ray echo 
should be detectable at a lower flux than an X-ray cluster, because of its small angular extent). 
The faintest X-ray flux detectable with the Advanced X-ray Astrophysics Facility should be of 
order 10“^® erg cm“^ s“^. Other future missions (such as XMM and Constellation X) have much 
larger apertures and could reach fainter fluxes, but with lower spatial resolution, so the photon 
noise is then dominated by the X-ray background for this very low flux. 


3. Effects of Gravitational Tensing 

The X-ray echo we have discussed should be visible when an intervening galaxy containing 
large amounts of dust is close to the line of sight to a gamma-ray burst. The gravitational field of 
the intervening galaxy will very often also deflect the light-rays in a very significant way. In fact. 
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the typical deflection angle of gravitational lenses caused by ordinary galaxies is of the same order 
as the scattering angles of the X-ray echos we have discussed in the previous section. 


In the presence of gravitational leasing, equation © giving the time delay needs to be 
modified to 

(1 + Zd) Dd D 


ctd = 


Dds 




(13) 


where (j) is the azimuthal angle. The projected potential il: is determined by = 2k, where 
the Laplacian operator differentiates with respect to the angular coordinates, k = Ti/Ticrit is the 
convergence, S is the surface density of mass, and Scrit = Ds/{DdDds) is the critical 

surface density (e.g., Blandford &: Narayan 1992). 


If the mass of the galaxy itself were negligible, and the gravitational leasing were caused by 
a smooth mass distribution around the galaxy (due to a possible cluster of galaxies of which the 
intervening dusty galaxy could be part of, or any large-scale structure density fluctuations present 
along the line of sight), then we can approximate the time delay surface with a second-order 
expansion. The images of background sources are then stretched along the two orthogonal axes of 
the shear by factors {1 — k — 7 )”^ and (1 — k -|- 7 )”^, where 7 is the shear (determined by the 
second derivatives of the projected potential). The X-ray echo would also be stretched by the 
same factors, becoming an ellipse. For high magnification, the X-ray echo increases to a given 
angular size in a much shorter time than in the absence of lensing, therefore increasing its flux 
through the factor eo/td in equation ( 12 ). 


When the mass of the galaxy is important, the gamma-ray burst may be multiply imaged. 
In this case, the X-ray echo should follow a curious evolution that could in principle provide a 
detailed map of the potential of the lensing mass. Considering a typical case of a three-image lens, 
with the third image being near the core of a galaxy and generally demagnified by a large factor, 
one should first see the image of the gamma-ray burst in the minimum of the time delay surface. 
This image would be followed by a highly elongated, arc-like X-ray echo expanding around the 
image. The arc would then increase in size and progressively bend around the center of the galaxy, 
until its two tips would touch on the opposite side of the first image. The second image of the 
gamma-ray burst should then arrive from the saddle-point of the time delay surface. The X-ray 
echo should then break up into two rings, one expanding out (with a shape that would continue 
mapping the mass distribution around the galaxy), and the other collapsing toward the center to 
eventually disappear in the third image of the gamma-ray burst (at the maximum of the time 
delay surface), expected to be very faint for the usual centrally peaked density profiles of galaxies. 


Gravitationally lensed X-ray echos should generally be brighter than non-lensed ones, because 
of the shorter time it takes for the echo to cover a given solid angle. In addition, the required 
scattering angle is reduced, and this can only result in an increased surface brightness for the echo 
(this increase is probably small, because most grains are small enough to produce scattering angles 
much larger than a few arc seconds). The observed time delays can again be used to obtain the 






product DdDs/Dds, as in other gravitational lenses. The advantage here is that, if the evolution 
of the X-ray echo can be closely monitored, the shape of the leasing potential can be mapped, 
eliminating the usual modeling uncertainties in other gravitational lenses. Notice, however, that 
the well-known “/t-degeneracy” is not eliminated; if a uniform sheet of matter with surface density 
S = Ko ^crit is added to the lens, resulting in the addition of ko9‘^I2 to the projected potential, 
and the potential ijj is multiplied by 1 — /to, the new time delay surface differs from the initial one 
only by the constant 1 — kq, which can be absorbed into the factor DaDg/D^s- Fortunately, this 
is the only systematic uncertainty that remains if the mapping of the X-ray echo can yield an 
accurate measurement of the shape of the time delay surface. 


4. Discussion 

Gamma-ray bursts having a dusty galaxy near their line of sight should be followed by an 
X-ray echo of the X-ray afterglow resulting from small-angle X-ray scattering by the dust. We 
have discussed the expected properties of these echos, including their angular extent as a function 
of time, and their flux. The main difficulty in detecting these X-ray echos will be their very low 
flux. Given the intensity of the observed X-ray afterglow following a gamma-ray burst, and the 
identification of any dust-rich galaxy near the line of sight from optical and infrared observations, 
it should be possible to predict approximately the flux of and time when the echo should be 
observable, and to estimate if it could be detected. 

There are some cases where the flux of the X-ray echo of a gamma-ray burst afterglow could 
be substantially higher than our estimate in §2. One possibility is the presence of gravitational 
leasing, which we have discussed in §3. Another possible case where a particularly bright X-ray 
echo might be observed is when the scattering galaxy is at very low redshift. In this case, the 
dust could have a larger angular extent, and the time delay should be much shorter for a given 
angle. For example, if (1 -|- Zd)Dd = 100 Mpc, then the angular size is 0 ~ 6" when the time delay 
is two months, and for the fluence cq = 2 x 10“®ergcm“^, the flux of the echo is now close to 
10“^^ ergcm“^ s“^. It needs to be pointed out that, given the presence of intervening dust, the 
probability density of the redshift of the intervening galaxy is constant at low redshift, because the 
cross section of a galaxy to produce an X-ray echo is proportional to its angular size (as opposed 
to the cross section for gravitational leasing, which is constant at low redshift for lenses with an 
isothermal profile). 

X-ray echos may also be searched around X-ray quasars with a dusty galaxy near their line 
of sight. In this case, the fluence cq should be the integrated variable flux over some feature in 
the lightcurve of the quasar, which can then be seen reverberated in the scattered echo with a 
measured time delay. Of course, the detection of the X-ray echo of a constant source would also 
be of great interest, even without the measurement of the time delay; in this case, high angular 
resolution is needed to separate the X-ray echo from the point spread function of the quasar, and 
any X-ray emission from a possible cluster around the quasar or intervening galaxy. 
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